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We set up a new calculational framework for the Yang-Mills vacuum transition amplitude in 
the Schrodinger representation. After integrating out hard-mode contributions perturbatively, we 
perform a gauge invariant gradient expansion of the ensuing soft mode action which renders a 
subsequent saddle point expansion for the vacuum overlap manageable. The standard "squeezed" 
approximation for the vacuum wave functional then allows for an essentially analytical treatment 
of physical amplitudes. Moreover, it leads to the identification of dominant and gauge invariant 
. . . classes of gauge field orbits which play the role of gluonic infrared (IR) degrees of freedom. Those 

emerge as a rich variety of (mostly solitonic) solutions to the saddle point equations which are 
' characterized by a common relative gauge orientation of the underlying gluon fields. We discuss their 

, scale stability, guaranteed by a virial theorem, and other general properties including topological 

quantum numbers and action bounds. We then find important saddle point solutions explicitly and 
}^^' examine their physical impact. Some of them are related to tunneling solutions of the classical 

, Yang-Mills equation, i.e. to instantons and merons, while others appear to play unprecedented 

roles. A remarkable new class of IR degrees of freedom comprises vortex and knot solutions of 
Faddeev-Niemi type, potentially related to glueballs. 

I. INTRODUCTION 
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fT^ , The strong couplings among soft QCD gluons manifest themselves in a variety of complex long-distance phenomena. 
' Most of them are thoroughly entwined with the vacuum state, as illustrated by such prominent examples as quark 
confinement, spontaneous chiral symmetry breaking, vacuum tunneling processes, the ensuing 9 structure as well 
. as large gluon condensates. Despite the apparent diversity of these and other effects, however, the essence of the 
■ underlying dynamics is often expected to involve just a few soft gluonic modes. 

. The quest for these long-wavelength excitations began shortly after the inception of QCD and has inspired the 
' development of various kinds of vacuum models, based e.g. on glueball condensation Q], Gaussian stochastic processes 
! 01 gluonic domains and instanton [5| as well as meron jg] ensembles. Over the last decade, lattice simulations 

I increasingly assisted in the search for predominant infrared (IR) gluon fields, mostly by means of numerical "filtering" 
Q-i! Q and gauge-fixing Q techniques. These simulations are now beginning to generate quantitative insights into the 
1^ ■ role of instantons and their size distribution 9] , and into the classic confinement scenarios based on (gauge-projected 
abelian) monopole ^ ilQj or center vortex {11^ J2J condensation. 

Nevertheless, at present no mechanism involving soft vacuum gluons can be uniquely or systematically related 
to QCD, and many crucial questions regarding the underlying fields, their stability, gauge-independent physical 
^ , interpretation, mutual interactions, relations to other vacuum fields etc. remain unanswered. Analytical progress in 
' this realm has been particularly slowed by the inevitable gauge dependence of the generally rather complex classical 
gluon field configurations on which most of the existing proposals are based. Similarly, approaches which reformulate 
non- Abelian gauge theory in terms of gauge-independent loop variables 13] or resolve the gauge constraints explicitly 
(e.g. in Coulomb gauge are often technically too involved for direct practical applications. 

In the present paper, we circumvent such complications by developing an approach in which not the contributions 
of single gauge fields but rather those of gauge invariant classes are treated jointly. These classes gather contributions 
from dominant gauge field orbits to low-energy Yang-Mills amplitudes and thus represent collective gluonic IR degrees 
of freedom. Technically, they are the saddle points of a soft mode action for gauge invariant matrix fields and therefore 
also provide the principal input for a systematic saddle point expansion of soft Yang-Mills amplitudes. 

Manifest gauge invariance is maintained throughout all calculations by working in the Hamiltonian formulation 
of Yang-Mills theory in the "coordinate" Schrodinger representation and by making use of explicit gauge projection 
operators. The Schrodinger picture is adopted mainly because it restricts gauge transformations to a fixed reference 
time, thereby effectively decoupling them from the dynamical time evolution, and because it often renders the impact 
of topological gluon properties particularly transparent (even without recourse to the semi-classical approximation) 
p.5]. For the reasons already alluded to, we will focus on gluonic effects and work in pure gauge theory without 
quarks. 

The individual IR degrees of freedom, i.e. the solutions of the saddle point equations, turn out to comprise a 
diverse range of specific features. Additionally, they have several important properties in common, including stability 
against scale transformations (an indispensable prerequisite for the saddle point expansion which is ensured by a 
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virial theorem) and characteristic topological properties inherited from the gauge group. The topology will turn out 
to be particularly useful for establishing relations between specific saddle-point families and the instanton and meron 
solutions of the classical Yang-Mills equation. 

Our approach maintains explicitly traceable links between the soft collective fields and the underlying gluon fields. 
The resulting IR dynamics is at an intermediate level of complexity, somewhere inbetwcen the microscopic theory 
itself and effective theories (e.g. for Polyakov loops) which just share the symmetries of the fundamental dynamics 
while coupling parameters have to be fitted to experimental (or lattice) data. Our soft-mode Lagrangian, in contrast, 
follows uniquely from the adopted vacuum wave functional and combines a reasonable amount of transparency with 
accessibility to essentially analytical treatment. 

The paper is organized as follows: in Sec. ^ we recapitulate the definition of the vacuum overlap amplitude in 
the Schrodinger picture. We then implement a gauge-projected vacuum wave functional on the basis of the Gaussian 
approximation and rewrite the overlap in terms of a bare action which had previously emerged in a variational context. 
In Sec. mil we take advantage of known 1-loop results to integrate out the hard-mode contributions to the bare action 
perturbatively. By means of a controlled derivative expansion, the renormalized soft-mode action density is then 
transformed it into a local Lagrangian which lends itself to direct analytical treatment. In Sec. IIVI we build on these 
results by establishing the IR-sensitive saddle point expansion for the functional integral over the soft modes and by 
deriving the saddle point equations whose nontrivial solutions constitute the new gluonic IR degrees of freedom. 

Important generic properties of these IR variables are established in Sec. including their scale stability due to a 
virial theorem, three topological quantum numbers and a lower bound of Bogomol'nyi type on their action. In Sec. IVII 
several classes of the more symmetric and most important saddle point solutions are found explicitly. They comprise 
topological soliton solutions of hedgehog type, which are related to classical solutions of the Yang-Mills equation, 
and solutions which carry different types of topological information and seem to have no obvious counterparts in 
classical Yang-Mills theory. One of the most interesting solution classes consists of solitonic links, twisted links and 
knots. Those emerge from a generalization of Faddeev-Niemi theory which turns out to be embedded in our soft-mode 
Lagrangian. In Sec. IVIII we classify all hedgehog soliton solutions, find their most important representatives explicitly, 
and establish the role of the regular solutions as mainly summarizing contributions from instanton and meron gauge 
orbits to the vacuum overlap. In Sec. IVIIII finally, we collect our principal results, comment on evaluating the 
contributions of the gluonic IR degrees of freedom to relevant amplitudes, and suggest directions for future work. 



II. VACUUM OVERLAP AMPLITUDE 



The vacuum overlap amplitude of SU (N) Yang-Mills theory (without matter fields) in the Schrodinger " coordinate" 
representation [3| reads 



Z' :^ (0,i+|0,t_) 



(1) 



The vacuum wave functional (VWF) depends on half of the canonical variables, i.e. on the static gauge fields 
A (x). Its gauge invariance, like that of any other physical state and wave functional, is dictated by GauB' law. This 
crucial requirement can be imposed on a given functional by simply projecting out its gauge-singlet component, i.e. 
by integration over the (compact) gauge group [l^ . Starting from an approximate and therefore generally gauge 
dependent wave functional ^o, one then obtains the associated VWF 



A" 



(2) 



{dfj, is the invariant Haar measure of the gauge group, Q is the homotopy degree or winding number of the group 
element U'''^\ and 9 is the vacuum angle) for which Gaufi' law is manifest. The vacuum energy has been set to zero. 

After interchanging the order of integration over gauge fields and gauge group in Eq. J^l, it becomes obvious that 
a gauge group volume can be factored out of Z' , i.e. 



Z' 



DU+ / DU- / DAiIj^ 



A'- 



Z / DU^ 



(3) 



since the A- integral is gauge invariant. In fact, the group volume Z' /Z is left over when the two un-normalized gauge 
projectors in the matrix element Z' are multiplied into one. The integrand of the remaining integral over the gauge 
group is naturally rewritten as a Boltzmann factor, i.e. 



Z= / DUei^v{-^b[U]), 



(4) 
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which defines the 3-dimensional Euchdean bare action Ti, as a functional of the "relative" gauge orientation U = 
UZ^U+ only. Owing to the gauge invariance of the gluon field measure, Ff, is gauge invariant as well and takes the 
explicit form 



[U] 



In 



V-o A 



(5) 



This action describes the dynamical correlations which the gauge projection of the functional t/jo in Eq. Q has 
generated. Hence it would become trivial, i.e. [/-independent, if ip^ were gauge invariant by itself. More specifically, 
Ff, [[/] gathers all those contributions to Z whose approximate vacua aX t — ±00 differ by the relative gauge 
orientation U . The variable U thus represents the contributions of a specifically weighted ensemble of all gluon field 
orbits to the vacuum overlap and is gauge invariant by construction. 

To proceed in an analytically tractable fashion, we now adopt the standard Gaussian approximation 



(G) 



A 



exp 



d'x / d'yA1{x)G-^'^''{x-y)J^,{y) 



(6) 



for the unprojected VWF [53, which has the decisive advantage of allowing integrals over A to be done exactly. As 

(G) 

expected from a ground state wave functional, i/iq has no nodes. It describes a "squeezed" state, i.e. an oscillator- 
type extension of the unstable coherent gluon states and thus the simplest natural candidate for the vacuum 
functional. In fact, Eq. © turns into the exact ground state for U (1) gauge theory (up to color factors) if the 
"covariance" is taken to be the inverse of the static vector propagator. Several additional properties indicate 
that Gaussian VWFs with suitably adapted covariances capture crucial features of the Yang-Mills dynamics as well. 
Indeed, with an appropriate choice for (see below) the wave functional © becomes exact at high momenta and 
incorporates asymptotic freedom. Moreover, it is known from variational analyses that Gaussian VWFs generate a 
dynamical mass gap and possibly confinement Il9l l20j . (Mass generation and most other features of 2-1-1 dimensional 
compact photodynamics are also reproduced |2l|.) Additional support for the Gaussian approximation will emerge 
from our results below. 

After specializing the expression ij^l for the bare action to "00*^^ j the functional integral over the gluon fields becomes 
Gaussian and can readily be carried out. The result takes the form of a 3-dimensional, bilocal nonlinear sigma model 



[U] = ^ 



d^x 



d'yLf {x)D^' {x-y)L\{y). 



(7) 



(Above we have omitted a term of higher order in the small bare coupling gi, which vanishes at the saddle points in 
which we will be interested below.) The [/-dependence enters F;, both via the one- forms 



L, (f) = [/t {x)d^U(x) 



(8) 



i.e. the Lie-algebra valued, left-invariant Maurer-Cartan "currents" (with real components L°), and through higher- 
order corrections to the bilocal operator 



{G + G^ 



lb 1 



(9) 



where G'^ = G"^ {x ~ y) [/^ (x) T^U (x) ® U (y) T''W {y), T'' = A°/2 and G-^"'' = G^M'^''. The above reformulation 
of the Yang-Mills vacuum overlap on the basis of a gauge-invariant Gaussian VWF was employed in Ref. as the 
starting point for a variational approach 56] . Alternatively, it can be obtained from a saddle point evaluation of the 
functional integral 23] in Eq. © which becomes exact for the Gaussian VWF jH3] • 

Although the nonlinear sigma model Q is easier to handle than the original Yang-Mills theory, its exact non- 
perturbative treatment remains beyond analytical reach p^ . Nevertheless, the parametric enhancement of the action 
(O by the factor g^'^ suggests that a useful approximation may be obtained from a saddle point expansion of the 
functional integral (0J). In order to render this approximation practical, however, one has to deal with the nonlocality 
of the bare action which encumbers the identification and evaluation of the saddle points. We will show in the 
following section that this can be efficiently accomplished by combining a renormalization group evolution of the bare 
action (which removes the explicit UV modes) with a subsequent derivative expansion to transform the IR dynamics 
into an approximately local soft-mode action. 
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III. SOFT GRADIENT EXPANSION 



For the reasons outlined in the introduction, we are mainly interested in soft Yang-Mills amplitudes with external 
momenta \pi\ smaller than a typical hadronic scale jjL (the lowest glueball mass, for example). This restricted focus 
permits us to recast the bare action {Tj) into a form which only retains soft field modes explicitly and which can be 
systematically approximated by a local Lagrangian. The present section describes the derivation and some useful 
features of this soft-mode action. 

Since the action incorporates asymptotic freedom (for proper choices of G, see below), the bare coupling gf, 
is small at the large cutoff scale Afjy where the theory is originally defined. The hard modes of the U field with 



momenta 



> /i can therefore be integrated out of the functional integral (0} perturbatively, down to values of the 

infrared scale /i where the renormalized coupling g (/i) ceases to be much smaller than unity. In practice, this may be 
done for instance by Wilson's momentum-shell technique ^25.j . after factorizing U into contributions from high- and 
low-frequency modes. To one-loop order, the resulting renormalization of the action just amounts to the replacement 
of the bare coupling gt by the running coupling g (/i) . This was confirmed in Ref. |2fil | where the one-loop coupling 
was obtained as 



' [N = 9b + In ■ 



(10) 



(for G (fc) = fc~^ at fc > /i). The scaling behavior of g (/i) makes asymptotic freedom explicit and reaffirms that the 
Gaussian VWF reproduces the qualitative UV behavior of Yang-Mills theory. In fact, Eq. IjlUI) equals the one- loop 
Yang-Mills coupling up to a small correction factor 1/11 which jirises from the absence of transverse gluons and could 
be avoided by introducing an anisotropic component for |2^. The one- loop integration over the high-momentum 
modes was found to be reliable down to /i ~ 1.3 GeV |2(i| which provides a useful benchmark for numerical estimates. 
In the valid range of /x values Eq. Q turns into the renormalized soft-mode action 



4g2 (^) 



d^x I d'yL'^^^,{x)G-'{x-y)L''^Jy) 



(11) 



where the subscript " < " indicates that U contains only k < fi modes. 

The action (|ll|l is still nonlocal. However, this nonlocality is substantially weaker than in the bare action (0 since 
the soft U fields in the integrand vary too little to resolve details of G^^ over distances smaller than /i^^. This 
observation can be turned into a controlled, local approximation scheme for the soft-mode action (|ll|l by exploiting 
the fact that the gradients of ?7< are bounded by the IR gluon mass scale. 



< /i- 



(12) 



Indeed, this bound suggests to expand the nonlocality of G ^ into derivatives which will act upon [/< after partial 
integration. Using the isotropy of G~-^ = G~^6ij (as mentioned above, an anisotropic component could be allowed in 
principle and would lead to somewhat more general expressions), one has 



G-i(x-y) 




(x-y). 



(13) 



The dimensionless constants Ci encode the low-momentum behavior of G ^ and could, e.g., be determined variationally. 
For our present purposes, however, it will be sufficient to adopt the standard expression G~^ (ji^ = \Jk'^+^'^ which 

approximates the solution of Schwinger-Dyson equations and variational estimates [Tol l20l | and incorporates both 
asymptotic freedom and a dynamical mass gap. The corresponding Ci can be read off directly from the Fourier 
transform 

1 ^l'K2{^i\x-y\l ^^^^ 



G-\x-y)^J-dl + iiH^ 



y) 



27r2 



\x - y\ 



{K2 is a McDonald function 27]), i.e. cq — 1, ci — —1/2, C2 = —1/8 etc. 

The combination of the above results leads to the intended reformulation of the nonlocal dynamics lO. As antici- 
pated, the bilocal action density for the soft modes in Eq. Hll|l becomes a ("quasi"-) local Lagrangian £ (x) and the 
action takes the familiar form 



r [U<] = J d^xC (f) . 



(15) 
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The Lagrangian is an expansion into powers of ^diU^ and therefore belongs to the class of generalized nonlinear 
sigma models. When expressed in terms of the Cartan-Mamer currents i<.i, it reads [s^ 

C (S) = -^^^J^i^ {^<.» (2^) L<,r {S) + ^5,L<,, (f) a,i<,, [x) - ^52l<,, [x] d^L^,, (x) + ...| . (16) 

We have omitted total derivatives AC from the higher-order terms of C since they do not affect the field equations. 
Nevertheless, they may generate non-vanishing surface terms due to infinite-action configurations which are generally 
irrelevant for the saddle point expansion. To lowest order, 

A/: (x) = 8^2^*^ [L<,^ {x) L<.. (f)] + ...} . (17) 

The gradient expansion in Eq. (|16|l is controlled by increasing powers of the parametrically small [dU^/ y)^. For 
practical purposes it can therefore be reliably truncated, at an order which is determined by the desired accuracy of the 
approximation to the exact action. Below, we will be interested in specific field configurations (saddle points) which 
we are going to find explicitly. For those one may directly check a posteriori whether the full action is sufficiently well 
reproduced and, if not, include contributions from higher gradients. Based on such tests and a virial theorem (see 

Sec. IV All , we found the truncation at O {^{d/ y)^^ to yield a generally sufficient approximation (at the few percent 

level) to the full action (|11|) . 

The first term in the Lagrangian (|16|l has the standard form of a 3-dimensional nonlinear cr-model (or principal 
chiral model). The second one, with four derivatives acting on C/<'s, is reminiscent of a similar term in the Skyrme 
model jl^. However, the identity 

1 2 

diLjdiLj = - [Li,Lj] +diLjdjL^ (18) 

shows that the four-derivative contribution to Eq. (|16|l contains, besides the commutator or Skyrme term, a part 
which qualitatively alters the character of the Euler-Lagrange equations. While the commutator generates only 
second-order terms to the field equations (see below), the piece without equivalent in the Skyrme model leads to 
additional fourth-order terms which allow for new types of solutions. (Several families of topological soliton solutions 
from the Lagrangian (|16|l will nevertheless turn out to resemble static Skyrmions jl^.) 

We end this section by emphasizing that the construction of an analogous gradient expansion in terms of the original 
gauge fields A would require the (residual) gauge freedom to be completely fixed and thus give rise to all the associated 
conceptual and calculational complications. (Otherwise a "soft" gauge field could just be turned into a "hard" one 
by a suitable gauge transformation and would spoil the "convergence" of the derivative expansion.) The locality and 
structural simplicity of the soft-mode Lagrangian (|16|l can therefore be regarded as a benefit of reformulating the 
dynamics in terms of the gauge invariant U field variables. 



IV. SADDLE POINT EXPANSION 



Our next task will be to set up the saddle point (or, more specifically, steepest descent) expansion of the functional 
integral over the soft modes, 

J W<exp(-r[(7<]), (19) 

where the vacuum overlap Z serves as the prototype for similar integrals in other soft amplitudes. This expansion is 
based on the IR saddle point fields Ui (x), i.e. the local minima of the soft-mode action (|15|l . The search for these 
minima is simplified by the fact that all finite-action t/-field configurations, including most saddle points, fall into 
disjoint topological classes (cf. Sec. IV Bl and below). Since fields which carry different topological charges - for now 
summarily denoted by Q - cannot be continuously deformed into each other, the local variation of the action may be 
performed in each topological sector separately. This amounts to solving the saddle point equations 

^U< [X) rj^^fj':'^) 



at fixed Q. 
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To leading order, the saddle-point expansion for the vacuum overlap Z can then be assembled by summing the 



contributions from the solutions u\^^ of Eq. (jJOJ, i.e 



V. 



(Q) 



exp 



(-r[f7P]) 



(21) 



where the pre-exponential factors Fi are generated by zero-mode contributions which typically arise if continuous 
symmetries of the action are broken by the solutions. The sum over the saddle points, labeled by i, symbolically 
includes integrals with the appropriate measure when the saddle points come in continuous families. Functional 
integrals for more complex amplitudes, including the gluonic Green functions, receive contributions from the same 
saddle points and are obtained by differentiating Z with respect to suitably implemented sources. Higher-order 
corrections to the approximation (|21|l can be systematically calculated from the (nondegenerate) fluctuations around 
the solutions uf^^. The reliability of the leading-order approximation (|21() increases with the action values of the 
saddle points because F [[/] ^ 1 prevents the saddle point contributions from being rendered insignificant by the 
fluctuations of 0(1) around them. In our case, compliance with this criterion is reinforced by the overall factor 
(/i) ^ 1 (for /i > 1.3 — 1.5 GeV) in the Lagrangian 1)16(1 which parametrically enhances the action (601] . 
For the explicit solution of the saddle-point equation (|20|l . as well as for part of our general analysis below, it is 
practically necessary to adopt a parametrization of the SU (N) group elements U which allows to work directly with 
their N'^ — 1 independent degrees of freedom. We will use the exponential representation 



[/< (x) = exp 



</>(f)n'^ (x)- 
2i 



(22) 



for this purpose, where the coefficient vector of the Lie algebra generators is decomposed into its direction, specified 
by the vector field n° with n°'n°- = 1 (which parametrizes the coset SU (N) /H where H is the Cartan subgroup of 
the gauge group), and its length </>. For simplicity, we will also specialize our following discussion to iV = 2. The 
soft-mode Lagrangian 1)16(1 can then be rewritten in terms of the unit vector field n° and the spin-0 field as 



C (U) = C2d ((/>, n) -I- Cm i<P, n) 



(23) 



where, for the reasons discussed in the paragraph below Eq. ((17II . all terms of the gradient expansion with up to four 
derivatives on the U fields are retained. The two-derivative part C2d^ i-e. the standard nonlinear cr-model, becomes 



C 



2d 



22^2 (^) 

while the four-derivative contributions turn into 



{d,(j)f + 2 {1- cos (j)) {d^h^f 



(24) 



C. 



1 



id 



2V (a*) A* 
-I- 4 sin (hdj 66 



d^(t)- sin (f>{d^h''y 
n^d^h" + 2(1- cos 



+ 2(l + cos(/)) id^(f>d^h'') 
(e''''"n''52r)% 2(l-c( 



2 (1 - cos(/>) {eijkd. 



2" 



(25) 



The general expressions above show that £ > 0, as required for the existence of the functional integral ((1911 . The 
same remains true if higher orders of the derivative expansion ((16(1 are included (as long as 5C/<//i < 1, of course). 
For the analysis of some generic saddle point solution properties, including their behavior under scale transformations 
(cf. Sec. 0), it will prove useful that both C2d and £4^ are even individually nonnegative. 

By varying the action ((15() of the Lagrangian ((23(1 with respect to (j) and n"', the saddle point equation l(2U() turns 
into a nonlinear system of four coupled partial differential equations. For one directly obtains 



^F [0, h] 
S(j3 (x) 



d-^(j) - 2(9^0 {dih^y - Ad^cjidAn^dkn" - 4 cos (/) {did jctydih" + djctyd^h" + difj^didju") dju" 



2 sin ( 



id,(t)din''y - (1 - COS0) {d,h''djn''y - 2d,n''d,d^n'' - [d^d^n'^y + [d^fi^y {d^n'') 



L cos {dS") {djfi'') - sin {S"'' + fi'^n'') d^h'^d^h'' - 2^^ q2^_ {d.h") 



= 0. 



(26) 



The three equations for the components of the n field, on the other hand, have to be derived by a constrained variation 
whose task it is to preserve the unit length of fi. As a consequence, they can be cast into the succinct form 



(27) 
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where the projection operator ensures that the action is only affected by those variations Sh" which maintain orthog- 
onality to The evaluation of the functional derivative with respect to h yields 



i5n° (x) 



sm( 



COS q)Oi<p 

(1 — COS 

2 



d'-n" ~2cos<j)^^<j)^J<j)^Jh'' 
dk<l}dkn'' + siiv(l)dk(j)d^dkn'' + sin^a.^ [d'^h" - h^ffd^n") 



sin (pdin!^ [din! 



n 

\2 



2 (1 — cosc 



(28) 



The saddle point equations (|26() and H27|l are independent of the coupling g because it enters the action only through 
the overall factor g~^. In general, their solutions have to be found numerically. However, we will demonstrate below 
that several nontrivial analytical solutions exist and that further important solution classes with a rather high degree 
of symmetry can be obtained by solving substantially simplified field equations. 

Moreover, essential qualitative solution properties can be derived without solving the saddle point equations ex- 
plicitly (cf. Sec. 0. Each topological charge sector contains at least one action minimum, for example, i.e. one 
solution of Eqs. H2t)|l and H27|) . The requirement of finding and including all of them would render the saddle point 
expansion practically useless. Fortunately, however, this turns out to be unnecessary. Below we will establish lower 
bounds on the action of the saddle points which are monotonically increasing functions of their topological charges. 
Hence contributions from saddle points in high topological charge sectors can generally be ignored. 



V. GENERAL PROPERTIES OF THE SADDLE POINT SOLUTIONS 

Before actually solving the saddle-point equations H2t)|) and H27|) in Sees. IVII and [VIII it will be useful to obtain a 
few general insights into the topological structure and stability properties of the solutions. This is the objective of 
the present section. 



A. Virial theorem 



In order to analyze the scaling behavior of the extended saddle point solutions and to establish the underlying virial 
theorem, we define the scaled fields 

0A:=0(Af), hl:^n''{Xx) (29) 

for real X ^ and note that /i is the only mass scale in the field equations (|26|1 and H27|) . This immediately implies 
that the solutions of the saddle point equations with scaled parameter /i — )■ A^^/i can be obtained by rescaling the 
original solutions (</>,«.) to {4>\,fix^. 

In the following, however, we will keep fi fixed. The scale-transformed extended solutions then cease to solve the 
field equations, and this simple observation together with two basic properties of the Lagrangian (|23|l can be turned 
into a virial theorem. The first step towards its derivation consists in establishing the relation between the actions of 
scaled and unsealed fields (as long as they stay finite), which can be read off from the 2- and 4-derivative parts of the 
Lagrangian H23|) separately: 

r (A) r fiA] = (A) + Tm (A) = jT^d (1) + Ar4rf (1) . (30) 

The second relevant property of the action based on Eq. (|23|l is its strict positivity for extended, i.e. not translationally 
invariant field configurations (cf. Sec. lIVp . 

r2d(i), r4d(i) >o. (31) 

The remaining step is to specialize the fields under consideration to the saddle point solutions. Since those extrem- 
ize the action under arbitrary small variations - which of course include infinitesimal scale transformations - one 
immediately has dV (A) /d\\-y^^ = and consequently the virial theorem 



T2d (1) = Tid (1) . 



(32) 
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Eq. H32() clearly exhibits the crucial role of the four-derivative terms in guaranteeing the stability of the saddle 
point solutions: for = the remaining nonlinear a-model action would be minimized by sending A — > oo (cf. Eq. 
(130(1 1. i.e. by the scale collapse of the solutions which Derrick's theorem predicts l29j. Our truncation of the gradient 
expansion therefore turns out to be the minimal approximation which can support localized, 

stable soliton solutions 61]. Furthermore, 

= 2T2d (1) > (33) 

A=l 

implies that the scaling extrema are indeed action minima and that solutions with (1) = are points of inflection. 

Finally, it is worth emphasizing that the coexistence of terms with different numbers of derivatives in the Lagrangian 
(|16|l . and thus ultimately the virial theorem H32I) and the existence of stable solutions, is brought about by the 
mass scale n which reflects the short-wavelength quantum fluctuations integrated out in Sec. IIIII This situation is 
reminiscent of the classical instanton solutions to the Yang-Mills equation whose typical size scale must likewise be 
generated by quantum fluctuations. 

To summarize, we have established a virial theorem which ensures that the extended solutions of Eqs. (|26|l and 
((27|1 are stable against scale transformations. In our context this is an indispensable property since unstable solutions 
would prohibit a useful saddle point expansion. As a side benefit, the virial theorem also provides stringent tests for 
the numerical solutions of the saddle point equations. 



d^T (A) 
dA2 



B. Topology 

As a three-dimensional principal chiral model with stabilizing higher-derivative terms, the soft-mode Lagrangian 
(|23|l allows for topological soliton solutions. In the present section we discuss three topological invariants or charges 
which such solutions and more general continuous fields U may carry jsj]. 

The most fundamental topological classification arises from the fact that all [/-fields with a finite action F based 
on any truncation of the Lagrangian l|16|l have to approach the same constant at — > oo. As a higher-dimensional 
analog of the stereographic projection, this compactifies their domain to a three-sphere S^p. All C/'s with finite F 
therefore describe continuous maps from into the "topologically active" part of the group manifold. For SU (N) 
the latter is the trivially embedded subgroup SU (2) ^ Sq. The resulting maps S^p — > Sq fall into disjoint homotopy 
classes, the elements of the third homotopy group tts (5'^) = Z, which are characterized by a topological degree or 
charge 

Q [U] - ^ J d^xe,,ktr {U^ d,UU^ djUU^ dkU] . (34) 

(The integrand can be shown to be a total derivative, as expected for a topological "charge density" .) For finite-action 
fields Q € Z. In terms of the ((/), n) parametrization ((22|l for U S SU (2), Eq. ((34|l reduces to 



1 



24^2 



d^x (cos </> - 1) e^jke''^''d^(t)n''djh''dkn''. (35) 



Two additional topological quantum numbers of U are carried solely by its n-field component. The first is the 
homotopy degree of the maps n° from the space boundary S^^ (where \x\ oo) into the unit sphere S\ on which fi°- 
takes values. Continuous maps of this sort are classified by the elements of the homotopy group 1^2 (S^) — Z. (The 
same topology characterizes the magnetic charge of Wu-Yang monopoles [23 ■) explicit integral representation of 
their degree is 



1 

8^ 



da^eijke'""'n''djn''dkn'' (36) 



where the integral extends over the closed surface dR^ — 5^ at \x\ — > 00. As expected, qm is ^-independent. 

The third topological invariant owes its existence to the fact that all n fields of finite action are required to approach 
a constant unit vector at spacial infinity. As above, this requirement compactifies into and thereby turns fi (x) 
into continuous maps S^p — s- S*?. Such maps carry a Hopf charge qh which labels the elements of the homotopy 
group TTa (<S'2) = Z, i.e. the Hopf bundle 33] . An explicit integral representation for qh can be constructed by means 
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of the local isomorphism between the nonlinear O (3) and CP (1) fields which expresses ri" in terms of a complex 
2-component field x — (xii X2) with unit modulus x^X — XiXi + X2X2 = 1 as = x^'''"X- Then one has 

ix)] = 4^ / d'xe^.k ix^d.x) {d.x^dkx) ■ (37) 

(A local integral representation for qn directly in terms of the n-field does not exist.) In contrast to the topological 
charges Q and qm which are of winding-number type, the Hopf invariant qn is a linking number. The underlying 
topological structure enables and classifies the link and knot solutions to be encountered in Sec. IVI CI 

Finally, we recall that the U field topology - as summarized in the conserved topological quantum numbers Q, q^ 
and qu - characterizes not only the saddle point solutions but a much larger class of continuous field configurations 
with finite and in some cases even infinite (see below) action. 



C. Bogomol'nyi bound 

The distribution of the solutions to Eqs. H2t)|) and (|27|l over a denumerably infinite set of topological charge 
sectors allows for the existence of more saddle points than could be handled in practical applications of the expansion 
(|2I|I . Hence additional criteria are required to select the most relevant saddle points in a controlled fashion. Such 
criteria will be established below, in the form of action bounds which are monotonically increasing functions of the 
absolute topological charge values. These bounds imply that contributions from saddle points with higher topological 
quantum numbers to functional integrals are increasingly suppressed by the Boltzmann factor exp (—T) and can be 
systematically neglected. In the present section we establish the action bound for fields which carry finite values of 
Q. A similar bound for fields with nonvanishing Hopf charge will be obtained in Sec. IVI CI 

The lower bound on the action of the Lagrangian (|23|l for any field U with a well-defined homotopy degree H34|l can 
be derived from the Lie-algebra valued expression 

Mi.j = adiLj + bsijkLk (38) 

where a,b G R and Li are the components of the Maurer-Cartan one- form ||SJ|. With the help of the Maurer-Cartan 
identity 

d,L,-d,L, = [L,,L,] (39) 

one obtains for its square 

M.jMij = a^d.Ljd.Lj - 2abe,juL^LjLk + 2b^L,L, (40) 

which obeys the basic inequality 

tr{A'UjM,j} <0. (41) 

(Recall from Eq. |(SJ) that the are expanded into anti-hermitean generators.) 

After specializing the bound (|4f |l to 6 = 1/^/2 and a = ±1/ [V^n) [b^, multiplying by — /i/ (2.g^), integrating over 
X and using the integral representation for Q, one arrives at 



J d^xtr + S.d,L,d,L,^ > T^^Q P] . (42) 



The more stringent of these inequalities results from the lower (upper) sign on the right-hand side if Q > (Q < 0). 
Their left-hand side amounts to the action which is produced by the first two terms [63 of the Lagrangian HI()|) . Hence 
the expressions (|42|l for both signs combine into an inequality of Bogomol'nyi type, 

r[U]>^^\Q[U]\. (43) 
9 (W 

This is the desired lower bound on the action of any sufhciently smooth U field with well-defined degree Q. It is 
saturated by those fields which obey the Bogomol'nyi-type equation 



(44) 
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where the lower (upper) sign again refers to Q > (Q < 0). The equations (|44|l can be considered as analogs of 
the self-(anti)-duality equations of Yang-Mills theory and have the interesting consequence that the rescaled Maurer- 
Cartan currents Li := ^Li/ (2«/x) of minimal-action fields in any Q-sector become generators of the su (2) Lie algebra, 



Li,Lj — iSijkLk. 



(45) 



Translationally invariant solutions (cf. Sec. IVI All . for example, have Li = and therefore trivially saturate the bound 
in the Q — sector. It remains to be seen whether nontrivial solutions in sectors with larger |Q| exist as well [iH 

The large factor multiplying IQI in the inequality H43(l indicates that contributions from saddle points with higher 
IQI are strongly suppressed. In fact, they seem safely negligible in most amplitudes which receive nonvanishing 
contributions from the Q — sector. However, one should keep in mind that even contributions with extremely small 
"fugacities" can sometimes have an important qualitative impact on the partition function. A case in point are the 
decisive monopolc contributions in the 2+1 dimensional Yang-Mills-Higgs model 35]. The physical interpretation of 
the IQI 7^ solutions and their analogs in Yang-Mills theory will be discussed in Sec. IVIII 



VI. IMPORTANT SADDLE POINT SOLUTION CLASSES 



In the following sections we are going to solve the four saddle point equations H26|) and H27() explicitly. As already 
mentioned, our main focus will be on solutions with a relatively high amount of symmetry since their typically smaller 
action values enhance their contributions to the saddle point expansion. Besides playing a predominant role in most 
amplitudes, these solutions can often be obtained either analytically or with moderate numerical effort. 



A. Translationally invariant solutions 



The simplest solutions of the saddle point equations and lf77|l are the x-independent matrices 



Uc 



exp 



2i 



const. 



(46) 



where cjjc and n° are both constant. These solutions form the complete vacuum manifold of the dynamics (|23|l . i.e. 
the set of all fields which attain the absolute action minimum 



r [Uc] = 0. (47) 

Due to a redundancy in the parametrization (|22() , the subset of vacua in the center of the gauge group is completely 
n°-independent : 

0e,fe = 2fc7r, C/e,fe - (-1)'^ . (48) 

In addition, those are the only vacua which do not break the global U (2) symmetry of the Lagrangian H23|l sponta- 
neously. A glance at the integral representation shows that none of the Uc carries topological charges Q ^ 0. 



B. Constant-n solutions 



Any constant vector n" solves the saddle point equation H27() identically and reduces the other one, Eq. (|26|l for 0, 
to the linear field equation 

{d^(f> - 2^i^(b) = 0. (49) 

Obviously, the solutions of this equation constitute families of new saddle points which differ by an additive constant 
and have degenerate action values. (This trivially ensures periodicity in (j>.) Alternatively, Eq. (|49|1 and the action of 
its solutions can be derived from the reduced Lagrangian [65j 

(50) 



£(n=c) 



1 



23r;2 (,,) II. 
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The presence of the 4th-order term in Eq. H49(l aUows for solution types which have no equivalent in Skyrme models. 
(Indeed, the commutator or Skyrme term (cf. Eq. (|18|l ') alone leads to a Laplace equation for 0, without a mass 
scale and with only constant regular finite-action solutions on i?'^.) A glance at the integrals (|35|l - (|37|l shows that 
solutions with constant h do not carry any topological quantum numbers, i.e. 

Q("-)=g(r^)=gg-=)=0. (51) 

All solutions of the linear field equation H49|) can be constructed by standard Green function techniques. The 
perhaps most straightforward approach is to fold the static Klein-Gordon propagator 



A(x-y;m) = - / — 3-.^^ = - — ^^-^ (52) 



with a "scalar potential" <& which is defined both to be a solution of the Laplace equation, 9^$ = 0, and to act as the 
inhomogeneity of the static Klein-Gordon equation 

a^c/. - 2^20 = (53) 

The field equation H49|l is recovered from Eq. (|53|l by applying the Laplacian to both sides. Since the potential $ plays 
the role of a source for the (p field, inversion of the Klein-Gordon operator immediately yields the general solution 

(x) = J d^j/A - y; V2^l^ $ (y) . (54) 

Of course, the regular finite-action solutions form but a small subset of those comprised in Eq. H54() . 

Spherically symmetric solutions can be obtained more directly by restricting the angular dependence of 0, i.e. by 
substituting the ansatz (/) (r) with r := \x\ into the general equation H49|) and ignoring for the moment potential 
singularities at the origin. This yields the radial equation 

r(/)"" + 4</)"' - 2^2 (^^/. ^ 2(1)') = (55) 
(radial derivatives d/dr are denoted by a prime) whose four linearly independent solutions can be found analytically: 



The subset of regular finite-action solutions is therefore of the form 



The associated potential 



(r) = ci + (l - e-^'^'- ) . (57) 

V2/i7' ^ 



shows that the expression (|57|1 in fact solves a generalization of the homogeneous field equation (|49|1 . with an additional 
delta-function singularity at the origin. Eq. H57|) is therefore a solution of Eq. H49(l everywhere except at = and, 
strictly speaking, one of its Green functions. A representative of this solution class is drawn in Fig. 1. 

After insertion into Eq. ()15|l. based on the Lagrangian H23|l. and use of the virial theorem H32|l one finds the solutions 
(|57|l to have the action 



r 





27r^ /■ 




" f (A^) io 



dr \ 



This action is not subject to topological bounds and reaches the absolute minimum F = for C2 = where the 
constant-n solutions turn into the translationally invariant vacua of Sec. IVI Al (In contrast to the center elements 
(|48|l . however, the value of (f> remains unrestricted here.) Due to their partly very small action values, the constant-n 
solutions may have a strong impact on the saddle point expansion which should be explored in detail by studying 
their contributions to suitable amplitudes. 
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C. Faddeev-Niemi type knot solutions 



In addition to the translationally invariant saddle points of Sec. IVI Al there are other and less trivial solutions of 
the field equations H2t)|) and (|27|l with constant (f> fields. Among them, the most intriguing class has the general form 



bc2,k = (2fc + 1) ^, U,2.k {x) = {-if ir'^n" (x) , 



(60) 



which satisfies Eq. (|26|l identically and carries no topological charge Q. In fact, the nontrivial topology of Uc2 (as 
that of any other constant-0 field configuration) has to reside exclusively in the a?-dependence of its n field, whose 
dynamics is governed by the (fc-independent) equation 



~abc 



Adjfi^dih'^djh'^ 



92 an" 



0. 



(61) 



The field equation (|61|l follows from the general saddle point equation (|27|l by substituting 4)c2.k and simultaneously 
plays the role of a (static) continuity equation for the conserved O (3) current. Alternatively, it can be obtained by 
directly varying the reduced Lagrangian 



9^ (m) 



1 



which follows from Eq. (|23(l after specialization to 0c2./c and reproduces the action (|15(l for the Uc2,k- 
the Lagrangian (|62|l is a generalization of the static Skyrme- Faddeev-Niemi (SEN) model [^E 1^3 



£_(SFN) 



2A2 



e 



(62) 

Remarkably, 
(63) 



In contrast to the SEN model, however, which was postulated on the basis of qualitative symmetry and renormalization 
group arguments [36L l37l| . our Lagrangian H62|l follows uniquely from the Yang-Mills dynamics and the Gaussian 
approximation to the vacuum wave functional. All coefficients are therefore fixed in terms of the IR scale /i and the 
coupling g (/i), i.e. Eq. H62|) does not contain free parameters. 

Particular solutions of equation (|61|l are n° = const., which belong to the class of translationally invariant vacua (cf. 
Sec. IVI B|l ■ and n° = (except at r = 0) which is an example from the "hedgehog" solution family whose detailed 
discussion will be the subject of the following sections. The (f)c2 hedgehog has infinite action (since U2c,k develops 
a monopole-type singularity at if = 0, cf. Sec. IVII C|l and its Lagrangian reduces exactly to the Eaddeev-Niemi 
form H63|l . Eq. (|61|l does probably also have cylindrically symmetric vortex solutions which are analogs of the "baby 
Skyrmion" solutions 40] in similar models. 

The most interesting and many- faceted solution classes of the field equation l|61() . however, are expected to be twists, 
linked loops and knots made of closed fluxtubes. Indeed, an intriguing variety of such topological soliton solutions was 
found numerically for the SEN part H63|l of the Lagrangian (|62|l in Refs. (371 141| . These solutions generally lack axial 
symmetry and carry a finite Hopf charge l66j 9^^'' 7^ 0. Moreover, their number and complexity increases strongly 



with the value of 



(c2) 

1h 



As in the higher- IQj solution sectors discussed previously, a practically useful saddle point 

expansion thus requires an effective means for selecting the relevant contributions in a controlled fashion. 

As anticipated in Sec. IV CI such a means can be provided by establishing that the action r'^'^'=^' based on the 
Lagrangian l|62|) is bounded from below by a monotonically increasing function of {qnl- Actually, this just requires a 
straightforward adaptation of a known bound on the SEN action [s^ • One combines the obvious inequalities 



> 



ff (a*) 

2 



(a,n") 



1 



(fix {d^fi") 



1/2 



(f'x{e''^''h''d^n''djn''y 



1/2 



(64) 



(65) 



(the first one holds because the omitted term in the Lagrangian H62|) is manifestly non-negative; the second one is a 



consequence of (a 



by > for any real a, b) with the Sobolev-type inequality |38 
1 



,3/2 



< 



26y233/47r4 
and a simple inequality due to Ward |3^ . 



fix^J{e'''"'h^^^n''^Jh^ f ( fix {e''^''fi°-dih^djn 



> 2 (e"'"' 



(66) 



(67) 
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to end up with the bound 

r[<Pc2,n]>^—^^\qH[n]f'. (68) 

(29/233/8^2 ^ 337^ rpj^ig bound is rather rough and could probably be made more stringent by incorporating the 
3rd term of the Lagrangian H62|l and by refining the Sobolev bound (|66|l which is expected to remain valid with about 
half of the factor on its right-hand side • 

The probably most important lesson of the present section is a new physical interpretation for Faddeev-Niemi-type 
knot solutions. In our framework, they reemerge as gauge invariant IR degrees of freedom which represent gluon 
field ensembles with a nonvanishing "collective" Hopf charge in the vacuum overlap and other amplitudes. This new 
interpretation may actually put the tentative identification of knot solutions with glueballs, advocated as a natural 
generalization of the fluxtube picture for quark-antiquark mesons in Refs. (37.. .42) . on a more solid basis. Indeed, 
the original association of Refs. js^, E3 is obscured by the interpretation of the n field as a gauge-dependent local 
color direction. Our -h field, on the other hand, is manifestly gauge invariant Is^. Moreover, glueballs are anyhow 
natural candidates for gluonic IR degrees of freedom in the Q — sector, so that their (perhaps partial or indirect) 
appearance in the saddle point solution spectrum would not be unexpected. 

An additional advantage of our new framework for the knot dynamics is that it allows the investigation of potential 
relations between the solutions of Eq. (|61|l and specific Yang-Mills fields which may play important roles in the 
vacuum, including e.g. topologically nontrivial pure-gauge fields in a non-linear maximally Abelian gauge '4^ and 
center vortices [ij. If such relations exist, they could perhaps be qualitatively traced by analytical methods. A full 
quantitative survey of the knot solution sector, however, will require a devoted numerical effort [ssf. 



D. Hedgehog solutions 



The existence of Skyrmions |2q in nonlinear cr-models with higher-derivative interactions suggests that our field 
equations (|26|l and (|27(l have topological soliton solutions of "hedgehog" type. 



n"(f)=^^ (f ) = (r) 



(69) 



(x" = x/r, r = \x\), as well. Their SU (2) "grand spin" symmetry characterizes the invariance of the corresponding U 
fields under simultaneous spatial and internal rotations and implies a substantial simplification of their dynamics. (For 
larger gauge groups SU {N) with N > 2, the three components of x form the part of the n° field which parametrizes 
the trivially embedded SU (2) subgroup.) In the present section, we discuss general properties of the hedgehog fields 
and derive their reduced Lagrangian and saddle point equation. In the subsequent sections I V 1 El and I V ill we find the 
most important solution classes explicitly and determine their physical interpretation. 

The principal topological characteristic of the hedgehog configurations Hti9|l is their 713 {S^) winding number Q. For 
fields of the form (|69() , its integral representation (|35|l reduces to 



(hh) 



1 

2^ JO 
1 

2^ 



(/.C*'')' (cos0('''*)-l) d,- 
sin (00) — sin 0'^'*''' (0) 



(0) - (C50) 



(70) 



(As expected from a topological invariant, it depends only on the boundary values of the field.) In Sec. IV Bl 
we established that finite action fields carry integer values of Q, and Eq. H70() confirms this explicitly. Indeed, the 
parametrization 1)22(1 for hedgehog fields (|69|l implies that well-defined U fields necessitate the boundary condition 
(jy{hh) j-Q-j _ 2A;i7r and that finite-action fields must additionally satisfy 0^'*''-' (00) — 2k2n (see below) where fci^2 and 
thus Q = ki ~ k2 are integers. Nevertheless, it is instructive to consider the more general boundary conditions 



^C''^) (0) 



Q 



(hh) 



m 



(71) 



(n, m, integer) which admit infinite-action fields with half-integer winding numbers as well (for either m or n odd) 
[69| . We will show in Sec. IVIII that hedgehog solutions to the saddle point equations (|2SJ), (|77|l under the boundary 
conditions (|71|l . both with finite and infinite action, can indeed be found. The hedgehog solutions with Q = ±1 will 
be of particular importance since they probably dominate all Q contributions to the saddle point expansion. 
This follows from the bound H43|) and from Skyrme-model type arguments |28| | which suggest that the minimal-action 
solutions in the IQj = 1 sectors are hedgehogs. 
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In addition, all hedgehog fields carry one unit of a second topological quantum number, the monopole-type charge 
Qm G ("S*^) • This becomes explicit when evaluating the {<j> independent) integral representation H3()|l for with 

qi^'^):=q^[S:]^l. J da^^^l. (72) 

Obviously, gm'*'' is independent of the boundary conditions for (j) and therefore of Q. The field with qm'*'' = — 1 (the 
"anti-monopolc" ) is obtained by replacing Xi by — ii, which corresponds to U ^ W for fixed 0. Eq. (|37|l reveals, 
finally, that the Hopf charge of all hedgehog configurations vanishes. 

The dynamics of the hedgehog fields is governed by the soft-mode Lagrangian (|23ll . Since x is an identical solution 
of the general field equation (|27|l for n, Eq. H23|l can be directly specialized to = x". Hence the integration over 
angles becomes trivial and the hedgehog action turns into 

/.no 

drC^'''^^ (r) (73) 

where is a 0'^'''' ■'-dependent radial Lagrangian. After substituting the ansatz (|69ll into the full Lagrangian H23|l . 

dropping total derivatives, suppressing the superscript of ^C*'*) and again denoting radial derivatives d/dr by a prime, 
one arrives at the explicit expression 





-j 




Jo 



9^ ifJ-) fJ- 



;^(r0")' + (3 + MV) (00 +4^2 



C0S( 



(74) 



All terms in are nonnegative. This has the consequence that each of them must vanish individually at any 

absolute action minimum. The complete set of hedgehog vacua is therefore (^4>l!^k^ j n'^'^'^^^ = {2kTr,x) and forms a 
subset of the translationally invariant center elements H48() . As anticipated, any finite-action solution of the form H69() 
has to approach one of these constant minima when r ^ oo. The constant solutions 4'[!2^k ~ (^^ -f 1) tt, on the other 
hand, are maxima of the action. A representative of this type was already encountered in Sec. IVI CI 

The radial equation for 0^'*'') can be derived by inserting the ansatz H69|l into the general field equation 1)26(1 or, 
more directly, by varying the radial Lagrangian (|74() . Either way, the result is 

r'^cj)"" + irct)'" - 2 (2 + ^V^) cj)" - ifi^rc/)' + 4^^ sin = 0, (75) 

i.e. an ordinary nonlinear differential equation of fourth order and of Fuchsian type. The associated boundary value 
problem can be solved numerically with rather moderate computer resources. The exploration of the full solution 
space is aided by two discrete symmetries of Eq. (|75|l which imply that any solution cf) (r) gives rise to the additional 
solutions —4> (r) and (f> (r) -I- 2mr. The former is a consequence of F [C/] = F [U^ while the latter simply reflects the 
periodicity in the angular variable cf). 

Not surprisingly, the field equation H75|) comprises the Gribov equation '44'| . It consists of the terms proportional to 
which originate from the nonlinear-cr-model part of the Lagrangian (|74|l and dominate when /i becomes the largest 
scale and/or when the higher derivatives become small |7(i. The analogy between the nonlinear potential term cx sin 
and a one-dimensional pendulum in a gravitational field |7l|. often used to characterize the solution spectrum of the 
Gribov equation, therefore applies to Eq. (|75|) as well. The stable (unstable) equilibrium positions of the "pendulum" 
are cj) = tt (0 = 0), modulo a multiple of 27r which represents additional full turns. It will be shown in Sec. IVIII that 
this analogy suffices to understand the qualitative behavior of all numerical solutions. 



E. Analytical hedgehog solutions by series expansion 

The (0, h) parametrization H22I) of the U field implies that regular solutions 0*^'*'') of the radial hedgehog equation 
(|75|l approach a multiple of 2tt at the origin. Their small-r behavior can therefore be determined analytically, either 
by expanding the nonlinearity of Eq. (|75|) into powers of small deviations Scj) (r) from the constant action minima 
'^c'fc'^ = 2kTT or by expanding the r dependence of the full solution into a Frobcnius series. Similarly, finite-action 
solutions can be obtained for r ^ oo by asymptotically expanding around the hedgehog vacua. Inside their regions 
of validity, these expansions provide useful insights into the qualitative behavior of the hedgehog solutions as well as 
quantitative checks on the numerical solutions to be found in Sec. IVIII 
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We start by deriving the solutions of the hnearized hedgehog equation and the corresponding power series expansion 
around the origin. Inserting the ansatz 

(/) (r) = 2fc7r + 5(j) (r) + O ((50^) (76) 

into the radial saddle point equation (|75(l and retaining only terms up to first order in 5(f), one arrives at the fourth- 
order linear differential equation 

r^54)"" + Ar5(t)"' - 2 (2 + ^?r^) 54>" - 4fi^r6(j)' + 4^i^dcf> = (77) 

which can be solved analytically by standard techniques. The general solution is a superposition of four linearly 
independent base solutions Sipi, 



>(r)=^Q<50a^), (78) 



1=1 



whose real, dimensionless coefficients Ci remain undetermined and have to be specified by imposing initial or boundary 
conditions. The base solutions {S(f)i} are 



5c^3= {l + V2fir)^^, 5c^4= {l-V2^ir) (79) 



The requirements of regularity and uniqueness on the solutions at the origin dictate two of the boundary conditions. 
The first one, Scj) (0) = 0, implies (p (0) = 2fc7r and thus ensures uniqueness at r = while the second one, 6(j)" (0) = 0, 
is then imposed by the behavior of the base solutions ((75)) . Accordingly, the small-r behavior of the general regular 
solution is restricted to 



4> (r) 2mr + cifir + C2ti^r^ + O (/r^) (80) 

where the constants Ci^2 are linear combinations of the 01^3^4 whose values can e.g. be specified by providing initial 
data for (/>' (0) and 4>"' (0). All higher-order coefficients of the expansion are then fixed. 

Alternatively, one can obtain the solutions of the full, nonlinear saddle point equation H75|l towards r — > by 
analytical continuation into a Frobenius series. A somewhat tedious calculation yields 



2-337 



A3 



2 i , 



103 



30/ 2^'^'' 5! 



+ O (r^) (81) 



where the coefficients (f>i and 03 are again left to be determined by initial conditions. Even-order derivatives of (or 
equivalently the coefRcients 02fc) vanish at r = while those of odd order, 02A;+i, can be expressed in terms of 0i and 
03. A comparison between Eqs. H78|) and H81|) shows that the solutions of the exact radial equation start to differ 
from those of the linearized equation at O (r^) . Hence the series solution (I8III permits a more accurate check of the 
numerical solutions over a larger radial interval. 

Analogous expansions around the constant action minima 0^\'''' — 2fc7r exist asymptotically, i.e. towards r — s- cx), for 
all finite-action solutions. Infinite- act ion solutions of Eq. (|75|l . finally, can be linearized around the constant solutions 
4>^2k^ ~ (2^ + 1) TT which they approach at one or both ends of the radial domain. The resulting equation for (50 differs 
from Eq. 1)77(1 only in the sign of the 50 term. Its solutions are linear combinations of generalized hypergeometric 
functions. 



VII. NUMERICAL ANALYSIS AND PHYSICAL INTERPRETATION 



We now turn to the numerical solution of the hedgehog saddle point equation (|75|l . Due to the periodicity in 0, 
the considered range of boundary values can be limited without loss of generality to (0) € ]0, 27r]. Regularity at the 
origin then further specifies (0) = 2tt and imposes 0" (0) = (see Sec. IVI E|) . The value of the topological charge Q 
fixes a third boundary condition, 

0(c5o) = 2^(1-Q), (82) 

owing to Eq. (|70|) . Hence all regular hedgehog solutions in a given Q-sector can be found by just varying the value 
of a fourth boundary condition. In the following, we use the initial slope /3 := 0' (0) for this purpose. At the end of 
the section, we will also find irregular solutions with (0) = tt. 
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A. Instanton classes 



We begin our exploration of the hedgehog solution space by searching for the regular finite-action solutions of Eq. 
(|75|l which, as established in Sec. IV Bl carry integer values of Q. The numerical analysis shows (and the pendulum 
analogy in Sec. IVII CI will explain) that only one solution of this type exists in each Q sector. In the simplest case, 
(5 = 0, this is just the translationally invariant vacuum solution 0^'''*' = 27r. 

For Q = 1 we find the prototypical nontrivial hedgehog solution, depicted in Fig. 2. In order to clarify its physical 
interpretation, we note that it shares the tts (S*^) homotopy classification, encoded in the topological charge Q \U\ of 
the relative gauge orientation U = UZ^U+, with the Yang-Mills instanton Of course, both also share the saddle 
point property, as the instanton minimizes the classical Euclidean Yang-Mills action in the Q = I sector. In order to 
trace their association further, we inspect the relative gauge orientation Ui^ym of a Yang-Mills instanton with size p. 
It is of hedgehog form as well, and its x-dependence in (Euclidean) temporal gauge is known to be _46j 

h,YM (r) ^== (83) 



\/r^ + p' 

in the parametrization (|69|l . For a direct comparison with oxlt Q ^ 1 solution, we have included (pi^vM with p = 
as the dashed curve in Fig. 2 (and adapted it to our periodicity interval convention by adding 27r). The radial 
dependence of both can be seen to be surprisingly similar. This implies that the dominant contributions from all 
Q — 1 gauge field orbits to the vacuum overlap have a relative gauge orientation close to that of the Yang-Mills 
instanton and indicates that the Q = 1 hedgehog solution primarily summarizes contributions from the instanton 
orbit. (Of course, one would not expect exact agreement since our solutions contain scale-symmetry breaking quantum 
corrections [t^ and contributions from other Q = 1 gauge fields as well.) Accordingly, and generalizing the above 
findings to multi-instanton solutions, we will refer to the unique regular finite-action solution of Eq. (|75|l with integer 
Q as the "Q-instanton class" . 

Our hedgehog saddle point equation H75I) and its instanton class solutions derive from the gradient-expanded 
soft-mode Lagrangian l)16|l . It is instructive to compare this approach with a variational estimate of one- instanton 
contributions to the bare action ((3J) in Ref. [47| |. By approximately minimizing the bare action with one-parameter 
families of trial functions similar to the instanton profile (|83|l and using qualitative scaling properties, it was argued 
in Ref. ^47] that radiative corrections can stabilize the instanton size. Our exact saddle point solutions make the 
dynamical size stabilization manifest. We have already traced the underlying mechanism to the virial theorem l|32|) 
which is independent of most specific features of the soft-mode dynamics and thus overcomes the chronic infrared 
instability of dilute instanton gases 4Sj in a rather generic way. For p ~ 1.5 GeV, the size p ~ 2/i~^ j73|| of the 1- 
instanton class solution agrees inside errors with the results of instanton liquid model ^| and lattice simulations. 
It also assures that the two leading terms of the gradient expansion (|13|) yield a sufficiently accurate approximation 
to the instanton action (cf. the comments below Eq. (|17|l ). 

Our 1-instanton class profile function 0/ (r) is rather similar to the one found by approximately minimizing the bare 
action JSJ variationally |42|| . This indicates that the bulk of the instanton's physics and size distribution is generated by 
soft modes, as one would intuitively expect. Our approach therefore provides a well-adapted and efficient framework 
for the treatment of these and other vacuum fields. In contrast to variational approaches, furthermore, it allows 
to systematically find all saddle points exactly (including those which are not of hedgehog form). Already in the 
hedgehog sector, for example, we will find solutions with more complex and unprecedented shapes than Eq. I|83|l . 
Since there is little guidance for the choice of suitable trial functions in these and other cases, such solutions would 
be difficult to find variationally. 

According to Eq. (|7T|l . all monotonic hedgehog solutions with Q > (Q < 0) have negative (positive) slopes 
(3 = (j)' (0) at the origin. The anti-instanton class with Q — —1, in particular, results from changing the sign of the 
instanton boundary value, /3j = and can be obtained without further calculation: it simply results from the 

combined action of the two symmetry transformations — > — and (/> — > 0-|-47r on the instanton class solution. Hence 
the Q = ±1 instanton classes have degenerate action values, precisely as their Yang-Mills counterparts. 

Multi-instanton class solutions are characterized by an integer topological charge Q > 2. The modulus \Pqj\ of 
their (negative) initial slope grows monotonically with Q, i.e. 

\Pi,Q'\>\Pi,q\ for Q'>Q- (84) 

As in the 1-instanton case, the action-degenerate multi-antiinstanton classes with Q < —2 can be constructed by fiip- 
ping the sign of the corresponding multi-instanton classes and adding An. The initial slopes of the multi-(anti)instanton 
solutions are therefore related by 



Pl-Q = ~Pl,Q- 



(85) 
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The IQI > 2 hedgehog solutions correspond to special arrangements of the underlying Yang- Mills multi- 
(anti)instantons. In fact, the relative gauge orientation U — UZ^U+ of a multi-instanton configuration is of hedgehog 
type only if all individual (anti)instantons are centered at the origin. This raises the question whether \Q\ > 2 
Yang-Mills instanton solutions with separated individual positions are at least approximately represented by other 
solutions of the saddle point equations (|26|l and H27I) . Experience from Skyrme-type models, whose analogous |Q| > 2 
Skyrmion solutions are well approximated by rational 52] or harmonic 53] maps, might suggest that similar types of 
non-hedgehog field configurations approximate higher-Q solutions of Eqs. H2t)|) and (|27|l as well. 

The action of the 1-instanton class solution is large (F/ ^ 220/g^ at /i = 1 GeV), in analogy with the large action of 
typical Yang-Mills instantons, and its direct impact on the saddle point expansion is therefore small 01 . Moreover, the 
action bound (|43|) implies that higher-Q instanton classes should be irrelevant for most amplitudes, with potentially 
important exceptions as mentioned in Sec. IV CI Instanton "liquid" vacuum models (ILMs) are built on the same 
premise and suggest that physically far more relevant contributions originate instead from ensembles of instantons 
and anti- instantons with equal average densities. It would be important to determine whether contributions of this 
sort are approximately represented by nontrivial Q = solutions of Eqs. and l(?7|l as well. In any case, our above 
results imply that they cannot be of hedgehog type. 

B. Meron classes 

We now extend our search to hedgehog solutions with infinite action. Although their relevance for the saddle 
point expansion is not obvious, one might speculate that their infinite-action suppression could be overcome by 
some additional mechanism (see below). Our chief motivation for discussing them here, however, derives from their 
association with the infinite- act ion meron solutions @ of the classical (Euclidean) Yang-Mills equation. 

Hedgehog solutions with infinite action are far more the rule than the exception. In fact, all regular (0 (0) = 2tt) 
solutions of Eq. (|75|l with initial slopes (3 inbetween the discrete set of instanton-class values have infinite action 
since they approach one of the constant fields (f>M (oo) = (2A: -I- 1) tt towards spacial infinity. The latter carry a nonzero 
action density (cf. Eq. (|74|l ) and furthermore imply that the corresponding, asymptotic U fields Um — *■ oo) remain 
angle-dependent. Exactly the same behavior characterizes the relative gauge orientations U = UZ^U+ of Yang-Mills 
merons in temporal gauge, which are of hedgehog form as well. Furthermore, Eq. (|82() shows that solutions with 
(j) (oo) = {2k -I- 1) TT carry half-integer topological charge Q, again as the Yang-Mills merons. In analogy with the 
instanton-class solutions of the previous section, we will therefore call these solutions "2(5-meron classes" . 

The profile function (f>M (f) of a typical 1-meron class solution with Q = 1/2 is drawn in Fig. 3. A direct comparison 
with the corresponding profile of the Yang-Mills meron in temporal gauge, 

4>M,YM {r) = lim 0/,YM {r) = -vr6' (r) (86) 

(where 9 is the step function), is complicated by the fact that the classical meron is pointlike while our solutions 
incorporate quantum effects which break dilatation symmetry and stabilize their size at a finite value. Such effects 
are expected to smoothen the singularity of the Yang-Mills meron, too, and probably cause our solution (j)M (r) to 
become non-monotonic by overshooting in the transition region. We therefore draw the Yang-Mills meron profile in 
Fig. 3 (dashed curve) with the finite size p = 2/i of the instanton class solutions (and adapt it to our periodicity 
interval convention). 

Our nomenclature for multi- meron classes with \Q\ > 1/2 includes only solutions with half- integer topological charge 
because solutions with even "meron number" and correspondingly integer Q coincide with the Q-instanton classes. 
This is expected since the relative gauge orientation Um,ym of the underlying Yang-Mills multi-merons is of hedgehog 
type only if all individual merons sit on top of each other. Such configurations, when carrying integer overall values of 
Q, coalesce into the corresponding Yang- Mills instantons and those are represented by the Q-instanton class solutions 
in our framework. 

The behavior of the multi-meron class solutions with Q > 3/2 is qualitatively rather similar to that of the 1- 
meron class, although size and strength scales may differ substantially. As an example, Fig. 4 shows a typical 
3-meron class solution. An important general property of all solutions with half-integer Q is that they come in 
families of continuously varying sizes. As already alluded to, this is because their size depends on a second mass scale 
/3a/ = (fi'jyj (0) (in addition to /i) and because solutions for all values of Pm in the finite intervals 

Pi,Q > Pm,q+i/2 > Pi,Q+i (87) 

(where Q > 0, i.e. /3 < 0, and f3i^Q=o = are implied) can be found. As in the instanton sector, multi-anti-meron 
classes with negative Q are obtained from the positive-Q solutions by changing their sign and adding An. 
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The variable mass scale /?m in the meron sectors implies that there are infinitely more meron than instanton class 
solutions. This makes it tempting to speculate that the meron "entropy" contributions to the weight function of 
functional integrals over U might be able to overcome the infinite- act ion suppression. If so, it would shed new light 
on the physical interpretation not only of our solutions but also of the Yang-Mills merons themselves, whose potential 
role remains controversial. Furthermore, it would suggest a modified saddle point expansion in which action and 
entropy are minimized jointly. These issues deserve further investigation. 

C. Singular hedgehog solutions 

Above we have classified all regular hedgehog solutions, i.e. those which satisfy the initial condition (0) — 2tt. 
We are now going to examine the remaining solution classes of the radial field equation H75|l . Its members share the 
alternative initial condition ^(0) = tt, may display a rather complex spacial structure and are characterized by a 
monopole-type singularity at the origin, i.e. they solve Eq. (|75() everywhere except at af = 0. 

In order to understand the qualitative behavior of both regular and irregular hedgehog solution classes from a 
common perspective, it is useful to elaborate on the analogy between the hedgehog equation (|75|l and the pendulum 
equation which was mentioned in Sec. IVIDI According to this analogy, the instanton classes correspond to exactly 
Q full turns of the pendulum, where the sign of Q indicates the direction of the rotation. The pendulum mass starts 
in the unstable equilibrium position at time t = \nr = —oo with just enough initial speed (3 = 4'' (0) to finally end up 
there again for t = \n.r = +oo. This analogy implies, in particular, that there is exactly one regular hedgehog solution 
for each integer Q and that the constant (/> = 27r is the only regular solution with Q — 0. The meron class solutions 
start from the unstable equilibrium position as well. However, their initial velocity (3 is insufficient for completing all 
turns in full. The last turn remains uncompleted, i.e. the pendulum swings back, oscillates around and finally settles 
into the stable equilibrium position. Hence all meron class solutions have half-integer Q and are non-monotonic. 

For the irregular solutions, on the other hand, the pendulum starts at the stable equilibrium position = tt. 
When not provided with sufficient initial speed f3 to complete a full turn, it just performs damped oscillations around 
(j) = IT. The corresponding solution has Q = and is depicted in Fig. 5. When |/3| is sufficiently large, however, the 
pendulum can perform Q full turns before settling into the stable equilibrium position. As a consequence, all these 
solutions have integer Q ^ and infinite action (cf. Eq. H74(l '). Pursuing the analogy further, one would also expect 
irregular solutions which have an initial value </>' (0) exactly as needed to end up at the unstable equilibrium position 
when r — > oo (after possibly completing a number of full turns). Such configurations would carry a half- integer 
topologically charge Q and a finite action. For obvious reasons they turn out to be highly sensitive to variations of 
the initial condition, however, and therefore difficult to establish numerically. 

In contrast to the instanton and meron classes, the singular hedgehog solutions do not seem to have obvious analogs 
among the solutions of the classical Yang-Mills equation. As opposed to the instanton (meron) classes, furthermore, 
the irregular integer-Q (half-integer-Q) solutions have infinite (finite) action. This allows for nontrivial hedgehog 
solutions with (3 = 0, which are necessarily irregular at the origin. The physical interpretation of all irregular 
hedgehog solutions and their relevance for the saddle point expansion remain to be clarified. 

VIII. SUMMARY AND CONCLUSIONS 

The main results of this paper are a practicable saddle point expansion for the Yang-Mills vacuum overlap amplitude 
in terms of gauge invariant, local matrix fields and the identification of new gluonic IR degrees of freedom in this 
framework. After adopting a gauge-projected Gaussian approximation to the vacuum wave functional, the IR degrees 
of freedom can be obtained explicitly as the saddle points of a soft-mode action which gather contributions from 
dominant gluon field families to soft Yang-Mills amplitudes and thus represent collective properties of the Yang-Mills 
dynamics. Since their gauge invariant definition makes no reference to specific amplitudes, furthermore, the IR degrees 
of freedom are universal. They provide both new structural insights into the organization of the low-energy Yang-Mills 
dynamics and the principal input for a systematic saddle point expansion of soft amplitudes. 

Our survey of the saddle point solution space uncovered a diverse spectrum of IR degrees of freedom which carry 
several topological charges with associated Bogomol'nyi-type action bounds and obey a virial theorem which guaran- 
tees their scale stability. Solutions with a relatively high degree of symmetry were obtained either analytically or with 
modest numerical effort. Since solutions of this type are generally characterized by small action values and hence play 
a dominant role in the saddle point expansion, we have investigated their properties in some detail. Besides trans- 
lationally invariant vacua and analytical solutions with a fixed relative gauge orientation, we have found topological 
solitons of hedgehog, (vortex) link and knot types. 
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Some of the IR degrees of freedom have a transparent physical interpretation directly in terms of the underlying 
gluon fields. The contributions from the gauge orbits of Yang-Mills instanton and merons, in particular, are gathered 
by saddle point fields of hedgehog type which share their (integer or half-integer) topological charge and represent 
vacuum tunneling processes in the Hamiltonian formulation of non-Abelian gauge theory. Although our saddle point 
solutions contain quantum effects and potentially relevant contributions from other gauge fields, those in the instanton 
class turn out to reproduce the relative gauge orientation between the in- and out-vacua of the Yang-Mills instanton 
rather closely. The finite extent of our meron solution classes, on the other hand, is generated by quantum effects 
which smoothen the singularities of the classical, pointlike Yang-Mills merons. Nevertheless, our meron classes turn 
out to share the infinite action of their Yang- Mills counterparts. 

Among all those IR fields which carry one unit of topological (instanton) charge, the single (anti-) instanton classes 
are expected to attain the minimal action value. As a consequence of the action bound, they should therefore dominate 
the saddle point expansion in all topological charge sectors. Similar configurations emerged in a variational treatment 
along with qualitative arguments in favor of their size stabilization. In our approach, this stabilization is manifest 
in the exact instanton class solutions themselves. In fact, their size turns out to be fixed at about twice the inverse 
IR gluon mass scale and agrees inside errors with instanton liquid model and lattice results. The underlying virial 
theorem and the soft gluon mass generation therefore provide new insight into the mechanism by which the chronic 
infrared diseases of dilute Yang-Mills instanton gases are overcome. 

The sizes of the meron class solutions with half-integer topological charge turn out to be of a more complex origin. 
Besides the dynamical gluon mass they depend on a second, variable mass scale which is encoded in a boundary 
condition. Hence meron classes exist within large and continuous size ranges and consequently form a far more 
extensive solution family than the instanton classes. This opens up the hypothetical possibility for their entropy to 
overcome their infinite action suppression in functional integrals. Such a mechanism would not only help to clarify 
the physical impact of our solutions but also shed new light on the still controversial role of the Yang-Mills merons 
themselves. In addition to the instanton and meron classes, finally, there exists a third class of hedgehog solutions 
which contains a monopole-type singularity at the origin. These irregular solutions can carry half-integer and integer 
(including zero) topological charges as well, and generally have infinite action. Hence their potential physical relevance 
seems to depend on the existence of additional mechanisms which could both smoothen their singularity and overcome 
their infinite-action suppression. 

Several other remarkable families of IR degrees of freedom turn out to be represented by topological solitons of 
Faddeev-Niemi type, i.e. by (potentially twisted) links and knots. In fact, our saddle point dynamics contains a 
specific generalization of the Faddeev-Niemi Lagrangian and shows explicitly how it is embedded in the Gaussian 
approximation to the Yang-Mills vacuum wave functional. This puts Faddeev-Niemi theory into a new perspective, 
as the effective dynamics of dominant sets of gauge field orbits with a collective Hopf charge, and provides the 
underlying unit-vector field with a manifestly gauge invariant meaning. The latter would well accord with the tentative 
interpretation of knot solutions as glueballs by Faddeev, Niemi and coworkers. This and other interpretations could 
be tested in our framework by directly evaluating the impact of the knot saddle points on suitable amplitudes, e.g. 
on glueball correlation functions. 

More generally, the saddle point expansion allows the systematic calculation of contributions from all relevant IR 
degrees of freedom to functional integrals which represent soft Yang-Mills amplitudes. First calculations of this type, 
focusing on fundamental vacuum properties including gluon condensates and the topological susceptibility, are under- 
way. In addition, our approach makes it possible to analyze the gauge field content of any IR variable individually 
by applying standard functional techniques to the integrals over gluon fields with which they are associated. Inves- 
tigations of this sort would not only provide further structural insight into specific IR degrees of freedom and their 
physical role but may also shed new light on the dynamical mechanisms by which soft gauge fields organize themselves 
into collective degrees of freedom. 

The diverse topological properties of the IR variables demonstrate that the gauge-projected Gaussian wave func- 
tional not only captures the homotopy structure of the gauge group but also implements "derivative" topologies 
which further characterize the saddle point families. Due to the typical robustness of such topological properties, 
the related results are expected to remain at least qualitatively valid beyond the Gaussian approximation. Moreover, 
the saddle point expansion engenders the means to test this expectation quantitatively, by mapping out limitations 
of the underlying vacTium wave functional in comparison with lattice data. Extensions of our framework to suit- 
able supersymmetric gauge theories would even permit analytical tests of this sort, e.g. by tracing vestiges of the 
monopole-based confinement mechanism in the vacuum functional. The insights gained from such investigations may 
also provide specific guidance for the development of improved collective-mode actions and consequently generate 
systematic corrections to the IR variables. 

Our IR saddle point expansion can be extended in several directions. A first important task would be a more 
exhaustive survey of the saddle-point solution space which should encompass potentially relevant approximate solu- 
tions. The extension to QCD proper requires the generalization to the gauge group SU (3) and the implementation of 
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quarks into the Gaussian wave functional, both of which pose no conceptual problems. Most topological properties, 
in particular, reside in the trivially embedded SU (2) subgroup of the full color group and will remain unchanged. 
A sufficiently complete treatment of the quark-gluon interactions and their impact on the effective action, however, 
appears to be more challenging. 

Our approach opens up a variety of directions for future research. Besides those already mentioned, it would for 
example be interesting to explore relations between the gauge-invariant IR degrees of freedom and gauge-dependent 
gluonic structures (monopoles, vortices, branes etc.) and amplitudes (e.g. the 2-dimensional nonlocal gluon condensate 
and Green functions) which appear in gauge-fixed formulations. Another useful endeavor would be the calculation of 
those higher-dimensional vacuum condensates which provide the principal input for the operator product expansion of 
glueball correlators. Duality sum rules could then link the contributions from different IR degrees of freedom with the 
low-lying glueball spectrum [s^, e.g. as a precursor and complement to a direct saddle-point evaluation of glueball 
correlation functions. 
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Figure captions: 



Fig. 


1: 


An example of the n = const, solution class. As in all following 
Other values of can immediately be accommodated by scaling 
interval of ^ restrict its initial value to </> (0) G ]0, 27r]. 


figures, we display the solution at n = 1 GeV. 
the r-axis. Our conventions for the periodicity 


Fig. 


2: 


The canonical 1-instanton class solution with Q = 1. 




Fig. 


3: 


A typical mcron-class solution with Q = 1/2. 




Fig. 


4: 


A 3-mcron class solution with Q = 3/2. 




Fig. 


5: 


An example for a nontrivial Q = hedgehog solution with monopole-type behavior at the origin. 
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